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1  J.  ABSTMAC  T 


Radar  Thomson  scatter  observations  of  the  midday  ionosphere  over  Randle  Cliff 
have  been  made  during  March-May  1971.  The  correlation  between  various  ionospheric 
parameters  has  been  determined.  Emphasis  has  been  placed  on  the  effect  that  magnetic 
activity  has  on  the  shape  of  the  electron-density  distribution.  In  general,  following  a 
magnetic  impulse,  the  F2  maximum  height  is  initially  raised  in  proportion  to  the  amount 
of  magnetic  activity,  is  subsequently  loweied,  and  eventually  returns  to  the  equilibrium 
position.  Also  the  F2  maylmum  density  and  the  total  electron  content  (below  600  km) 
docronso  as  the  length  of  time  between  the  magnetic-activity  impvilso  and  the  measurement 
incrunscs.  This  situation  exists  for  at  least  a  day;  thereafter  the  diminution  is  reduced, 
and  the  profile  gradually  assumes  its  avoraco  shape. 
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ABSTRACT 


Radar  Thomson  scatter  observations  of  the  midday  ionosphere 
over  Randle  Cliff  have  been  made  during  March-May  1971.  The  corre¬ 
lation  between  various  ionospheric  parameters  has  been  determined. 
Emphasis  has  been  placed  on  the  effect  that  magnetic  activity  has 
on  the  shape  of  the  electron-density  distribution.  In  general,  follow¬ 
ing  a  magnetic  impulse,  the  F2  maximum  height  is  initially  raised  in 
proportion  to  the  amount  of  magnetic  activity,  is  subsequently 
lowered,  and  eventually  returns  to  the  equilibrium  position.  Also  the 
F2  maximum  density  and  the  total  electron  content  (below  600  km) 
decrease  as  the  lengto  of  time  between  the  magnetic-activity  impulse 
and  the  measurement  increases.  This  situation  exists  for  at  least  a 
day;  thereafter  the  diminution  is  reduced,  and  the  profile  gradually 
assumes  its  average  shape. 


PROBLEM  STATUS 

A  final  report  on  one  phase  of  the  problem.  Work  continues  on 
other  phases. 
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ELECTRON-DENSITY  PROFILES  OF  THE  IONOSPHERE  OBSERVED  NEAR 
WASHU^GTON,  D.C.,  DURDIG  THE  SPRING  OF  1971 


INTRODUCTION 

This  report  describes  the  results  of  radar  Thomson  scatter  observations  made  in  March, 
April,  and  May  1971.  An  interim  report  outlining  the  March  results  has  been  published  (1). 
The  ionospheric  response  to  magnetic  activity  was  of  particular  interest,  and  this  topic  will 
receive  considerably  more  attention  in  this  report.  The  experimental  configuration  and 
the  data-analysis  techniques  employed  in  the  determination  of  the  electron-density  pro¬ 
files  presented  here  have  been  reported  earlier  (1,2)  and  will  not  be  repeated. 

The  principal  results  discussed  in  this  report  have  been  presented  at  the  Fall  meeting 
of  the  American  Geophysical  Union,  and  all  of  the  profiles  have  been  made  available  to 
interested  scientists  through  the  facilities  of  the  World  Data  Center  (3). 


MEASUREMENTS 

A  significant  modification  has  been  made  in  the  method  of  coniputing  electron-density 
profiles  since  the  publication  of  the  March  results.  Instead  of  assuming  a  constant  value 
for  the  M  factor  (M  =  H  cos  6  sec  x).  its  full  altitude  dependence  has  been  incorporated 
into  the  computer  program.  This  modification  eliminates  a  systematic  error  of  as  much  as 
15%  in  the  profiles. 

Also,  each  of  the  nineteen  March  profiles  were  separately  fitted  to  a  fourth-degree 
least-squares  polynomial  to  deduce  the  various  ionospheric  parameters.  In  this  report  pro¬ 
files  were  obtained  by  combining  the  data  corresponding  to  the  same  day  of  observation 
prior  to  determination  of  the  polynomial.  In  addition  the  peak  F-region  density  (Np2max) 
for  each  day  has  been  taken  to  be  the  figure  dictated  by  the  polynomial  rather  than  the 
maximum  value  associated  with  the  set  of  unsmoothed  raw  estimates  of  density.  This  ap¬ 
pears  to  bring  the  value  of  Np2max  niore  in  line  with  the  Wallops  Island  ionosonde  results. 


Composite  Profiles  for  March,  April,  and  May 

Figure  1  is  a  scatter  plot  of  data  obtained  during  March,  April,  and  May  1971.  A 
least-squares  fourth-degree-polynomial  fit  to  the  raw  data  points  is  plotted.  The  decrease 
in  Np2max  between  March  and  May  is  seen  quite  clearly.  This  would  be  expected  if  past 
measurements  of  total  content  are  used  as  a  gauge,  for  it  is  well  known  that  the  total  con¬ 
tent  is  greater  at  the  equinox  than  during  the  summer.  Note  that  the  data  points  in  Fig. 

1  appear  to  have  preferred  locations  in  some  instances.  This  is  not  a  real  effect;  rather  it 
is  a  result  of  the  analysis  procedure  which  allows  only  quantum  jumps  in  the  Faraday  ro¬ 
tational  derivative.  This  has  been  discussed  in  an  earlier  report  (1). 
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profiles  each  are  based  on  eight  days  of  operation. 
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Table  1  lists  the  days  and  times  corresponding  to  the  individual  profiles  on  which  the 
curves  in  Fig.  1  are  based.  The  March  profile  is  based  on  ten  days  of  observation,  and  the 
April  and  May  profiles  are  based  on  eight  days  of  observation.  The  daily  profiles  are  given 
in  Appendix  A. 


Table  1 

Days  and  Times  for  the  Profiles  of  Fig.  1 


March  1971 

April  1971 

May  1971 

Day 

Time  (EST) 

Day 

Time  (EST) 

Q 

7'ime  (EST) 

12 

1346 

21 

1418 

11 

1450 

12 

1408 

21 

1444 

11 

1500 

15 

1450 

22 

1406 

■■ 

1353 

15 

1502 

22 

1418 

|H 

1353 

16 

1337 

23 

1411 

13 

1418 

1C 

1403 

23 

1425 

13 

1450 

17 

1310 

26 

1408 

14 

1440 

17 

1321 

26 

1418 

14 

1452 

18 

1334 

27 

1314 

17 

1451 

18 

1415 

27 

1345 

17 

1515 

19 

1320 

28 

1402 

19 

1338 

19 

1335 

28 

1434 

19 

1400 

22 

1320 

29 

1355 

20 

1502 

22 

1347 

29 

1438 

20 

1516 

23 

1607 

30 

1334 

21 

1330 

30 

1345 

21 

1339 

24 

1539 

24 

1551 

25 

1503 

25 

1518 

Ionospheric  Parameters 

Parameters  which  have  been  deduced  include  the  F2-maximum  density  (Np2max)» 
altitude  of  the  F2  maximum  (hp2max)»  electron  content  to  an  altitude  of  600  km 
[C(600)l ,  and  the  slab  thickness  to  an  altitude  of  COO  km  [t(600)]  .  These  parameters  are 
given  in  Table  2  for  each  day  of  operation,  and  the  average  values  and  standard  deviations  are 
given  in  Table  3.  The  upper  limit  for  the  computation  of  both  .slab  thickness  and  content  is 
dictated  by  the  rapid  diminution  of  Faraday  rotation  with  height  above  the  F2  maximum. 
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The  practical  upper  limit  was  ^600  km;  this  value  was  used  for  convenience  and  consist¬ 
ency*  as  a  termination  height  in  the  analysis. 


Table  2 

Atmospheric  Parameters 


Date 

(1971) 

NF2niax 

(10^^  electrons/m^) 

^F2max 

(km) 

C(600) 

(10^ electrons/m2) 

7(600) 

(km) 

March  12 

15 

310 

290 

3.23 

2.45 

232 

210 

16 

10,92 

280 

2.41 

221 

17 

12.00 

310 

2.82 

235 

18 

8.78 

300 

2.17 

248 

19 

13,76 

320 

3.23 

235 

22 

13.04 

290 

2.77 

212 

23 

13.67 

300 

2.85 

208 

24 

9.23 

310 

2.28 

247 

25 

8.22 

290 

1.83 

222 

April  21 

11.81 

340 

2.99 

254 

22 

7.80 

290 

1.80 

230 

23 

4.57 

310 

1.35 

296 

26 

9.25 

310 

2.26 

245 

27 

9.06 

290 

1.95 

216 

28 

6.30 

300 

1.48 

235 

29 

7.16 

300 

1.83 

256 

30 

7.15 

280 

1.71 

239 

May  11 

5.23 

300 

1.39 

265 

12 

8.39 

310 

2.26 

280 

13 

4.86 

320 

1.47 

302 

14 

5.85 

350 

1.80 

308 

17 

5.62  i 

360 

1.18 

210 

19 

5.44 

300 

1.43 

264 

20 

6.89 

300 

1.80 

262 

21 

7.04 

1 

280 

1.85 

263 

Table  3 

Average  Values  (A)  and  Standard  Deviations  (S.D.)  of  the  Ionospheric  Parameters 


Month 

and 

Year 

Peak  Electron 
Density 
(NF2max) 

F2  Maximum 
Height 
(hF2max) 

Electron  Content 
to  600  km 
IC(600)] 

Slab  Thickness 
to  600  km 

[T(600)] 

March 

1971 

A  =  11.52  X  10“ /m3 
S.D.  =  2.51xlO“/m3 

A  =  300  km 

S.D.  =  2.47  km 

A  =  2.60x  10“/m^ 
S.D.  =  0.46  X  10“/m2 

A  =  227  km 

S.D.  =14.64  km 

April 

1971 

A  =  7.88x  10“ /m^ 
S.D.  =  2.17xl0“/m3 

A  =  302  km 

S.D.  =  12.80  km 

A  =  1.92  X  10“/m2 
S.D.  =  0-.51xl0l7/m2 

A  =  246  km 

S.D.  =  23.87  km 

May 

1971 

A  =  6.15x  10“ /m3 
S.D.  =  1.18  xl0“/m3 

A  =  315  km 

S.D.  =  27.25  km 

A  =  1.65x  10“ /m2 
S.D.  =  0.34xl0“/m2 

A  =  269  km 

S.D.  =  30.1  km 
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The  average  times  corresponding  to  the  Randle  Cliff  Radar  (RCR)  observations  in 
March,  April,  and  May  were  1418, 1407,  and  J427  EST  respectively.  Table  4  shows  the 
average  values  of  hF2max»  f^2max.  and  the  total  content  below  600  km  [C(600)]  for  the 
Randle  Cliff  Radar  (RCR)  site.  For  comparision,  estimates  of  hF2max  'md  NF2max  have 
kindly  been  supplied  by  Dr.  J.  Nisbet*  of  Pennsylvania  State  University  on  the  basis  of 
the  CCIRt  model.  Also  shown  are  the  Wallops  Island  ionosonde  values  of  NF2max»  which 
are  based  on  7-E  form  data  supplied  by  Mr.  R.  GrayJ.  The  RCR  values  of  NF2max  are 
somewhat  higher  than  Nisbet’s  values  for  March  and  May  but  are  lower  in  April.  In  addi¬ 
tion  the  RCR  values  of  hF2max  are  consistently  higher  by  >20  km  than  those  obtained  in 
the  model.  The  Wallops  Island  values  of  NF2max  are  in  close  agreement  with  the  RCft  ob¬ 
servations,  except  during  May  when  they  differ  by  «10%.  The  percentage  deviation  is  less 
than  1%  in  March  and  2.8%  in  April. 


Table  4 

Average  Value  of  Ionospheric  Parameters  for  the  RCR  Site 


Parameter 

Month  and  Year 

Data  Source 

March  1971 

April  1971 

May  1971 

hF2max(km) 

300 

302 

315 

RCR  Thomson  Scatter 

273 

281 

281 

Penn  State  Model 

NF2max(m"3) 

11.52  X  lOll 

7.88  X  10“ 

6.15  X  1011 

RCR  Thomson  Scatter 

8.31 

8.22 

5.50 

Penn  State  Model 

11.49 

7.66 

5.50 

Wallops  Island* 

Total  Content 

2.60  X  1017 

1.92  X  1017 

1.65  X  1017 

RCR  Thomson  Scatter 

(m-2) 

2.29 

2.48 

1.78 

Penn  State  Model 

♦Wallops  Island  data  are  based  on  hourly  7E-fortn  data  at  1400  EST.  These  data  are  said  to  be  more  accurate 
than  the  F-plot  data. 


Crosscorrelation  of  the  Ionospheric  Parameters 

Table  5  shows  how  the  parameters  NF2max>  hF2max.  C(600),  and  t(600)  are  correlated 
during  March,  Apnl,  and  May  respectively.  The  behavior  during  March  and  April  is  quite 
similar,  b  t  the  May  behavior  is  considerably  different.  During  March  and  April,  fluctua¬ 
tions  in  the  F2  maximum  height  are  positively  correlated  with  NF2maxi  C(600),  and  t(600), 
but  these  fluctuations  are  negatively  correlated  with  the  same  p,irameters  during  May. 
Though  the  correlation  betv/een  NF2max  and  t(600)  is  negative  daring  the  first  two  months, 
there  is  no  correlation  between  those  two  parameters  during  May.  Also  there  is  little  corre¬ 
lation  between  C(600)  and  r{600)  during  March  and  April,  but  during  May  the  corre^tio’ 
is  significant  (+0.41).  The  only  consistent  correlation  is  that  linking  the  total  content  ind 


♦Private  communication. 

tCosnite  Consultatif  International  doi  Radiocommunications. 
iPrivate  communication. 
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the  F2  maximum  density.  The  correlation  was  observed  to  be  positive  in  all  three  cases, 
being  +C.85  in  March,  +0.84  in  April,  and  +0.77  in  May.  Thus,  as  expected,  NF2max 
strongly  controls  the  total  content.  The  exact  for  the  difference  betv/een  the  com¬ 
bined  March- April  behavior  and  that  for  May  is  not  known  at  present,  but  it  is  felt  that 
the  folio  win?/  facts  may  be  significant: 

•  The  peak  electron  density  and  total  content  for  May  were  lower  than  for  the 
two  previous  months.. 

•  The  monthly  averagfj  3-hourly  Fredericksburg  K  index,  KpR,  over  a  24-hr  period 
preceding  the  median  observation  time  was  fo>md  to  be  1.80  in  March,  2.12  in 
April,  and  2.05  in  May.  However,  K  indices  are  roughly  logarithmically  related 
to  magnetic  activity.  Since  a  minor  geomagnetic  storm  occurred  in  May  (and 
not  in  March  or  April),  one  finds  that  the  actual  magnetic  activity,  in  terms  of 
r*,  is  greater  in  May  than  in  either  of  the  other  two  months. 


Table  5 

CroivS'correlatinn  of  the  Ionospheric  Parameters 


Month/ParametJt-r 

t(600) 

C(600) 

hF2max 

March  ' 

-0.37 

+0.85 

+0.29 

April 

’■  ^''F'max 

-0.37 

+0.84 

+0.44 

May 

-0.01 

+0.77 

-0.27 

March 

1 

+0.54 

+0.51 

April 

■  hF2max 

+0.40 

+0.60 

May 

-0.11 

-0.29 

March 

-0.07 

April 

rC(600) 

-0.14 

May 

+0.41 

MODIFICATION  OF  IONOSPHERIC  PARAMETERS  BY  MAGNETIC  ACTIVITY 

To  determine  the  relationships,  if  any,  between  magnetic  activity  and  the  parameters 
NF2max»  hF2max»  t(600),  and  C(600),  a  crosscorrelation  function  was  constructed  between 
these  parameters  and  KpR.  Tables  6,  7,  and  8  list  the  lag  intervals,  the  times  to  which  they 
correspond,  and  the  associated  K  indicesf  for  March,  April,  and  May  respectively.  The 
equivalent  r  values  are  also  given. t 

*The  parameter  r  is  a  measure  of  the  disturbance  magnitude.  At  a  magnetic  latitude  of  50  deg  (close  to  the 
RCR  latitude)  the  association  to  be  made  between  K  and  r  is  approximately  logarithmic,  i.e,,  K  a  log  r. 
(See  K.  Davies,  Ionospheric  Radio  Propagation,  National  Bureau  of  Standards,  Monograph  80,  GPO,  Wash¬ 
ington,  D.C.,1965,  pp.  26-25.) 

tThe  listed  values  are  actually  sums  of  three  adjacent  KpR  indices  (which  are  basically  3-hourly).  These  new 
indices  will  te  symbolized  by  SKpR.  Each  I  KpR  index  thus  covers  a  9-hr  period  of  time,  and  every  third 
index  is  completely  independent. 

tThe  r  values  are  based  on  the  effective  three-hourly  average  KpR.  That  is,  one  finds  the  value  of  r  corre¬ 
sponding  to  one  third  the  tabulated  value  of  SKpR.  The  following  concspondencc  is  used: 

r:  0  5  10  20  40  70  120  200  330  500 

2Kfr/3:  012345  6  7  8  9 

and  linear  interpretation  is  used  for  values  of  £Kfr/3  which  fall  between  the  tabular  values. 
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Functions  which  describe  the  time-varying  correlation  between  the  3-hourly  Kfr  in¬ 
dices  and  the  parameters  NF2max.hF2max»  C(600),  and  t(600)  are  given  in  Figs.  2-4  cor¬ 
responding  to  March,  April,  and  May  respectively.  Prior  to  processing,  the  Kfr  indices 
were  smoothed  by  forming  a  running  avers^e  of  three  consecutive  three-hourly  indices. 

Thus  each  value  of  Krr  used  in  the  analysis,  and  tabulated  in  Tables  6-8,  is  actually  char¬ 
acteristic  of  a  9-hr  interval  of  time  and  is  denoted  by  the  sy  ’'ol  SKfr.  Nevertheless  the 
spacing  between  the  listings  in  Tables  6-8  is  only  3  hr,  since  the  basic  Kfr  data  is  3-hourly. 
The  lag  associated  with  each  9-hr  interval  is  actually  the  temporal  midpoint  of  that  inter¬ 
val;  thus,  for  example,  the  so-called  fifth  lag  interval,  which  extends  between  the  lag  times 
of  13.25  hr  and  22.25  hr,  is  assigi^ed  a  mean  lag  time  of  17.75  hr.  In  addition  the  corre¬ 
lation  coefficient  between  SKfr  and  an  ionospheric  parameter  for  the  fifth  lag  interval 
represents  the  ionospheric  response  to  an  impulse  of  magnetic  activity  which  occurred 
17.75  hr  earlier  on  the  average. 


Fig.  2— Correlation  functions  relating  the  dependence  of  the  iono¬ 
spheric  parameters  on  magnetic  activity  for  March  1971 


Fig.  3— Correlation  functions  relating  the  dependence  of  the  iono¬ 
spheric  parameters  on  magnetic  activity  for  April  1971 
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Fig.  4— Correlation  functions  relating  the  dependence  of  the  iono¬ 
spheric  parameters  on  magnetic  activity  for  May  1971 


Assuming  that  a  positive  excursion  in  magnetic  activity  (+5EKfr)  has  occurred,  then 
a  positive  correlation  coefficient  linking  SKfr  and  an  ionospheric  parameter  P  implies  that 
the  parameter  will  experience  a  positive  excursion  6P  >  0,  where  6P  is  proportional  to  the 
magnitude  of  the  correlation  coefficient.  If  the  correlation  were  negative,  then  a  positive 
excursion  (+5 SKfr)  implies  a  down  ward  excursion  in  the  parameter  5P  <  0,  where  again 
5P  is  proportional  to  the  magnitude  of  the  correlation  coefficient. 

Fig.  5  was  prepared  to  compare  more  readily  on  a  monthly  basis  the  effect  of  magne¬ 
tic  activity  on  a  particular  parameter,  say  NF2max*  The  monthly  comparisons  for  C(600), 
t(600),  and  hF2max  are  given  in  Figs.  6-8. 


Fig.  5— Correlation  functions  relating  the  dependence  of  NF2max 
magnetic  activity  for  March,  April,  and  May  1971 
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LAG  INTERVAL 


Fig.  6— Correlation  functions  relating  the  dependence  of  C(600)  on 
magnetic  activity  for  March,  April,  and  May  1971 


LAG  INTERVAL 


Fig.  7— Correlation  functions  relating  the  dependence  of  t(600)  on 
magnetic  activity  for  March,  April,  and  May  1971 

On  inspection  of  Fig.  5,  it  is  rather  obvious  that  the  F2  maximum  electron  density  is 
negatively  correlated  with  the  magnetic  activity.  Generally  speaking  the  correlation  at 
small  time  lags  is  low  but  slightly  positive  during  March  and  April  and  somewhat  negative 
during  May.  At  later  times  the  correlation  is  again  low  during  March  and  May  but  strongly 
negative  during  April.  One  consistent  feature  is  that  the  correlation  is  rather  strongly  nega¬ 
tive  between  lag  intervals  of  3  and  6  or  for  an  average  lag  of  «16.25  hr.  Hence  the  phe¬ 
nomena  which  are  responsible  for  the  diminution  of  electron  population  near  the  F2  maxi¬ 
mum  are  maximized  at  roughly  16  hr  following  an  impulse  of  magnetic  activity. 

Since  NF2max  and  C(600)  are  highly  correlated,  it  is  not  surprising  that  the  same  state¬ 
ments  just  made  concerning  NF2max  and  2Kfr  also  hold  the  C(600)  and  2Kfr.  One  can 
notice  the  similarity  of  Figs.  5  and  6. 
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Fig.  8— Correlation  functions  relating  the  dependence  of  hF2max 
magnetic  activity  for  March,  April,  and  May  1971 


As  seen  in  Fig.  8,  the  F2  maximum  height  hF2max  is  generally  positively  correlated 
with  SKfr  for  small  lag  times.  One  can  recall  that  the  actual  monthly  magnetic  activity, 
as  indicated  by  the  parameter  r,  satisifes  the  inequality  rMay  >  ^Aprii  >  rMarch>  where  the 
bar  denotes  the  average.  For  the  first  two  1^  intervals,  the  degree  of  correlation  between 
hF2max  and  SKfr  (p(hF2max.  SKfr,  lag))  exhibits  the  same  monthly  ordering,  i.e.,  PMay 
>P  April  >PMarch‘  Again,  although  the  situation  is  not  easy  to  generalize,  the  overall  pat¬ 
tern  could  be  described  as  follows;  First,  the  F2  maximum  height  is  raised  in  proportion 
to  the  amount  of  magnetic  activity.  Second,  the  F2  maximum  descends  below  its  equilib¬ 
rium  value;  this  condition  occurs  near  the  fifth  lag  interval  (~i7.75  hr)  in  April  and  near 
the  eighth  lag  interval  (~26.75  hr)  in  March  and  May.  Finally,  there  is  a  tendency  for  the 
F2  maximum  to  return  to  its  equilibrium  position  after  a  long  time. 

Figure  7  shows  how  the  so-called  slab  thickness  depends  on  the  lag  time  for  March, 
April,  and  May.  In  April  and  May  the  initial  response  is  an  increasing  negative  correlation 
between  t  and  SKpR  with  a  tendency  toward  more  positive  correlation  at  later  times.  The 
situation  is  reversed  in  March. 


Translation  of  the  Coefficient  of  Correlation  to  the  Fractional  Change  in  the 
Ionospheric  Parameter 


The  sample  correlation  coefficient  p  between  two  sets  of  data  A  and  B  is  defined  by 


P(A,  B) 


N 


— ^(Ai-A)(B,-B) 


(1) 
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where  and  as  are  the  standard  deviations  of  A  and  B  and  where  A  and  B  are  the  aver¬ 
ages.  Also  if  A  and  B  are  normally  distributed  variables,  it  is  well  known  that  the  line  of 
regression  of  A  upon  B  is  given  by 

A(B)  =  p{A,B)^B  .  (2) 

Ur 

This  equation  says  that  given  a  particular  value  of  B,  the  average  valu^of  A  is  determined 
(4).  This  is  not  to  be  confused  with  the  sampled  population  average  A.  Since  the  value 
of  A  is  dependent  on  B  in  a  linear  way,  a  small  change  (5B)  in  B  can  be  related  to  a  re¬ 
sultant  change  (6A)  in  A.  Taking  B  to  represent  the  set  of  SKpR  values  and  A  to  repre¬ 
sent  the  ionospheric  parameter  P,  the  following  very  useful  expression  can  be  deduced: 

^  =  p(P,IKfr)^=-^ - SSKfr  .  (3) 

r  r  OSKpR 

It  will  now  be  of  interest  to  calculate  the  fractional  response  of  an  ionospheric  parameter 
to  an  impulse  in  the  3-hourly  index  Kfr  .  For  an  average  jump  of  unity  in  the  3-hourly 
Kfr  index,  the  recorded  cumulative  index  SKpR  must  increase  by  a  factor  of  three.  If 
the  3-hourly  index  were  to  jump  by  an  amount  SKfr  ®  4,  then  6(!2Kfr)  ~  12.  Table  9 
gives  the  fractional  change  in  the  ionospheric  parameters  NF2max>  hF2max»  and  t(600) 
which  occurs  as  the  result  of  a  jump  SSKpR  =  12.  This  corresponds  to  a  9-hr  average 
jump  in  KpR  of  4.  It  is  noteworthy  that  the  fractions  listed  in  the  table  are  only  indica¬ 
tive  of  fluctuations  which  would  be  induced  in  the  average  ionosphere  by  an  impulse  of 
magnetic  activity  at  the  time-lag  index  specified. 

ILLUSTRATION  OF  MAGNETIC-ACTIVITY  RESPONSE 

A  useful  computer  program  has  been  constructed  from  which  it  is  possible  to  visualize 
more  easily  the  effect  of  magnetic  activity  on  the  ionosphere.  The  output  of  the  program 
is  a  series  of  plots  corresponding  to  the  shape  of  the  ionospheric  electron-density  distribu¬ 
tion  for  specified  time  lags.  The  plots  are  not  exact;  they  are  idealized  Chapman-like  func¬ 
tions  corresponding  to  three  parameters:  NF2max.  hF2max>  and  the  scale  height  H.  The 
parameters  NF2max  and  hF2max  are  measured  directly,  but  H  is  deduced  from  the  relation 

T  =  4.133  H,  (4) 

which  is  strictly  valid  only  for  a  Chapman  distribution  of  the  form 


N  =  NF2max  exp 


(5) 


where  Z  =  (h  -  hF2max)/H.  For  each  month  of  observation  an  idealized  Chapman  distri¬ 
bution  is  computed  having  parameters  NF2max.  hF2max »  and  H.  The  perturbed  distribu¬ 
tion  parameters  are  given  by 
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Table  9 

Fractional  Change  in  Ionospheric  Parameters  for  a  Jump  in 
in  Magnetic  Activity  of  52Kfr  =  12* 


Lag 

(m) 

March 

April 

May 

NF2max 

hF2max 

H,t 

NF2max 

hF2max 

H,r 

NF2max 

hF2max 

H,r 

2.76 

+0.26 

+0.06 

+0.01 

+0.18 

+0.15 

+0.01 

-0.14 

+0.16 

-0,15 

5.75 

-0.06 

+0.01 

+0,02 

-0.11 

+0.68 

-0.04 

-0.12 

+0.14 

-0,13 

8.75 

-0.33 

-0.02 

+0.05 

-0.29 

-0.27 

-0.13 

-0.14 

+0.13 

-0.17 

11.75 

-0.34 

-0.01 

+0.07 

-0.55 

-0.13 

-0.08 

-0.18 

+0.14 

-O.IG 

14.75 

-0.29 

-0.03 

+0.04 

-0.43 

-0.15 

-0.11 

-0.22 

+0.16 

-0.16 

17.75 

-0.23 

-0.04 

+0.02 

-0.41 

-0.12 

-0.06 

-0.30 

+0.27 

-0.21 

20.75 

-0.35 

-0.07 

-0.02 

-0.40 

-0.12 

-0.05 

-0.45 

+0.26 

-0.57 

23,75 

-0.47 

-0.09 

-0.02 

-0.41 

-0.11 

-0.06 

-0,11 

-0,19 

-0.23 

26.75 

-0.46 

-0.10 

-0.04 

-0.29 

-0.07 

-0,02 

-0,05 

-0,17 

-0.04 

29.75 

-0.32 

-0.07 

-0.03 

-0.43 

-0.09 

-0.03 

-0.04 

-0.12 

-0.03 

32.75 

-0.15 

-0.06 

-0.04 

-0.69 

-0.06 

-0.13 

-0.03 

-0.11 

-0.03 

35.75 

-0.09 

-0.04 

-0.03 

-0.62 

-0.07 

-0.19 

-0.04 

-0.08 

-0.04 

38.75 

-0.09 

-0.04 

-0.02 

-0.49 

-0.05 

-0.19 

-0.07 

-0.07 

-0.05 

41.75 

-0.05 

-0.05 

-0.02 

-0.42 

-0.01 

-0.18 

-0.20 

-0.05 

-0.08 

44,75 

-0.07 

-0.06 

-0.03 

-0.53 

-0.01 

-0.21 

-0,15 

-0.10 

-0.13 

♦52Kfr  =  12  corresponds  to  a  3-hourly  jump  in  KpR  of  4;  this  3-hourly  value  persists  for  9  hr  or  for  three 
successive  3'hourly  periods. 


NF2max  =  NF2max 


f  SNF2max 
I  NF2max 


hF2max  ~  hF2max 


and 


WF2max  “ 


t(600) 


4.133 


+  1 


f  Shp2max  \  ^  ^ 
\  hF2max  / 


(6A) 

(6B) 


Sr(600)  \ 
r(600)  / 


(6C) 


where  the  fractions  within  the  parentheses  were  obtained  in  the  previous  section.  Figures 
9-11  show  idealized  distributions  for  March-May  1971.  The  solid  curve  in  each  figure  is 
the  average  monthly  distribution,  and  the  dashed  curves  represent  the  perturbed  distributions. 


FSg.  0 — Idealized  response  of  the  ionosphere  to  magnetic  activity  as  a  function  of  the  time  lag  corresponding  to  the  Thomson  scatter  observations.  March 
data  are  depicted.  A  value  of  SSKpR  =  18  was  assumed.  This  corresponds  to  a  severe  jump  in  KpR  of  6  over  a  9-hr  period. 
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DISCUSSION 

It  has  been  kn  >wn  for  some  time  that  magnetic  storms  produce  an  effect  on  the  con¬ 
centration  of  electrons  in  the  ionosphere.  The  general  behavior  is  one  in  which  the  varia¬ 
tion  in  electron  density  is  an  increase  followed  by  a  decrease.  The  negative  main  phase  of 
the  magnetic  storm  is  normally  associated  with  a  drop  in  the  total  electron  content.  It  has 
been  shown  by  Jacchia  (5)  and  Newton  et  al.  (6)  that  the  neutral  species  are  heated  by 
magnetic  activity.  It  is  anticipated  that  this  heating  and  the  related  turbulence  produced 
in  the  atmosphere  will  increase  the  height  at  which  diffusive  separation  appears,  thus  in¬ 
creasing  the  electron  loss  rate  (7).  On  the  other  hand,  Taylor  (8)  has  observed  large  de¬ 
creases  in  electron  production  during  the  negative  phase  of  magnetic  storms— a  fact  which 
suggests  that  a  decrease  in  the  rate  of  electron  production  may  be  at  least  as  important  as 
an  increase  in  the  electron  loss  rate.  The  first  report  of  a  large-scale  enhancement  in  the 
total  electron  content  during  the  initial  positive  phase  of  a  magnetic  storm  was  made  by 
{joodman  (9)  using  synchronous-satellite  data  obtained  in  1967.  The  effect  has  been  stud¬ 
ied  rather  exhaustively  by  Papagiannis  et  al.  (10)  who  suggest  that  positive  storm  effects 
are  primarily  a  dusk  phenomenon.  Goodman  (9)  has  suggested  that  electrodynamic  forces 
may  play  a  role  in  the  enhancement  of  total  content,  and  this  possibility  has  been  given 
some  support  by  Evans  (11).  Recently  Jones  (12)  and  Jones  and  Rishbeth  (13)  have  stud¬ 
ied  the  storm-time  variation  of  the  F2-layer  electron  concentration  and  the  possible  origins 
of  the  variation.  They  claim  that  the  two  effects  are  produced  by  competing  processes 
and  that  the  positive  effect  is  produced  by  storm-induced  changes  in  the  thermospheric 
wind  pattern.  They  also  suggest  that  an  increase  in  the  equatorward  neutral  wind  will  drive 
the  F2  layer  to  greater  heights,  thus  increasing  the  electron  concentration,  since  the  loss 
coefficient  is  an  exponentially  decreasing  function  of  hei^t.  (Of  course,  the  neutrals  only 
move  horizontally,  but  their  motion  imparts  an  effectively  upward  movement  to  the  ions 
because  of  the  presence  of  the  magnetic  field.) 

The  present  data  were  obtained  during  periods  of  generally  quiet  geomagnetic  activity. 
Nevertheless,  one  minor  geomagnetic  storm  did  occur  during  May.  The  effects  of  this 
storm  may  be  seen  by  comparing  Figs.  A19-A26  in  the  appendix.  The  primary  thrust  of 
this  study,  however,  is  directed  toward  the  effect  of  low-level  magnetic  activity  on  the 
ionosphere.  The  depletion  in  electron  concentration  following  a  moderate  excursion  of 
magnetic  activity  is  illustrated  in  Fig.  12.  An  approximate  30%  drop  in  the  F2  peak  elec¬ 
tron  density  occurred  on  April  22  following  an  increase  in  magnetic  activity  on  April  21. 

In  addition  the  F2  layer  height  decreased  substantially  between  the  two  days.  It  is  antic¬ 
ipated  that  the  phenomena  responsible  for  this  behavior  are  active  for  all  degrees  of  mag¬ 
netic  activity.  Indeed,  the  results  of  an  analysis  of  the  Thomson  scatter  data  obtained  at 
Randle  Cliff  are  consistent  with  the  notion  of  a  depletion  in  electron  concentration  at 
some  time  following  a  small  impulse  of  magnetic  activity.  This  electron-content  depletion 
is  controlled  by  a  depletion  of  the  F2  maximum  density,  and  this  effect  is  generally  great¬ 
est  whenever  the  magnetic  activity  precedes  the  profile  measurement  by  «16  hr,  on  the 
average.  Such  a  lag  is  also  characteristic  of  the  ionospheric  response  to  large  magnetic 
storms.  The  immediate  response  of  the  ionosphere  (the  response  for  small  lag  times)  is  a 
slight  enhancement  in  the  total  content  during  March  and  April  and  a  decay  in  May.  The 
March-April  behavior  is  characteristic  of  the  initial  buildup  of  electron  content  which  occurs 
during  the  p'>sitive  phase  of  many  magnetic  storms.  Even  during  May  the  amount  of  initial 
electron-content  decay  was  less  at  the  early  lag  times  than  during  later  times  (a:15  hr). 
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Evidently  two  competing  processes  are  in  progress,  with  the  positive-phase  process  occurring 
first  and  the  negative-phase  process  occurring  somewhat  later.  It  is  possible  that  the  nega¬ 
tive-phase  process  is  broader  during  May  and  adds  a  negative  bias  to  the  initial  positive 
response.  During  March  and  April,  the  negative  phases  are  not  as  broad. 


It  is  emphasized  that  the  results  previously  presented  in  the  section  on  the  modifica¬ 
tion  -of  parameters  assume  normally  distributed  raiidom  variables.  It  is  anticipated  that  the 
variables  (KpR  and  the  ionospheric  parameters  NF2ma.\,  hF2max.  C(600),  and  t(600))  will 
be  so  distributed  in  nature.  Unfortunately  the  populations  are  sampled  so  seldom  for  each 
lag  interval  (between  eight  and  ten  times)  that  it  is  not  possible  to  establish  this  fact  with 
any  statistical  significance.  Naturally,  by  combining  the  runs  from  all  three  months,  it 
should  be  possible  to  conduct  significant  chi-square  tests  to  determine  the  goodness  of  fit 
for  each  distribution.  Unfortunately  such  a  combination  of  the  three  months  is  not  appro¬ 
priate,  due  to  the  rather  obvious  seasonal  dependence  in  the  data.  A  possible  approach  to 
solving  this  problem  would  be  to  remove  the  seasonal  bias  prior  to  processing.  This 
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approach  has  not  been  attempted  as  of  this  writing.  It  will  be  carried  out  eventually,  ho«.'- 
ever,  and  the  results  will  be  contained  in  a  final  report. 


SUMMARY 

Radar  Thomson  scatter  observations  of  the  midday  ionosphere  were  made  during  March, 
April,  and  May  1971.  From  these  observations  the  following  F-region  parameters  were 
deduced:  the  F-region  peak  density  NF2max.  the  F2  maximum  height  hF2max.  the  total 
content  below  600  km  C(600),  and  the  equivalent  slab  thickness  t(600).  The  average 
monthly  electron  population  decreased  between  March  and  May,  as  expected,  since  the 
equinoctial  ionosphere  is  more  dense  than  the  summer  ionosphere  at  midlatitudes.  It  was 
found  that  the  total  content  is  strongly  controlled  by  the  peak  F-region  density  for  all 
three  months.  In  addition,  fluctuations  in  the  F2  maximum  height  were  positively  corre¬ 
lated  vrith  fluctuations  in  the  NF2max>  C(600),  and  t(600)  during  March  and  April;  the  con¬ 
verse  was  true  in  May.  During  May  there  was  almost  no  correlation  between  the  param¬ 
eters  Npamax  and  t(600),  however,  during  March  and  April  the  correlation  was  negative. 

It  was  also  found  that  C(600)  and  r(600)  were  almost  independent  during  the  first  two 
months,  but  during  May  they  were  positively  correlated.  It  is  suggested  that  the  differ¬ 
ences  in  the  March-April  and  the  May  results  may  be  partially  explained  by  the  differences 
in  the  character  of  the  magnetic  activity  during  the  two  periods. 

An  examination  of  the  time  lag  between  measurements  of  the  Fredericksburg  Kfr 
index  and  the  Thomson  scatter  observations  leads  to  the  construction  of  crosscorrelation 
functions  between  the  degree  of  magnetic  activity  and  Vfjrious  ionospheric  parameters.  As¬ 
suming  a  gaussian  dish-ibution  of  the  random  variables,  the  average  ionospheric  response  to 
an  impulse  of  magnetic  activity  has  been  deduced.  The  general  behavior  is  one  in  which 
the  F2  maximum  height  increases  initially  in  proportion  to  the  amount  of  magnetic  activ¬ 
ity,  subsequently  decreases,  and  eventually  returns  to  its  equilibrium  value.  It  is  found 
also,  that  the  F2  maximum  density  and  the  total  content  decrease  with  increasing  lag  time, 
the  greatest  diminution  occuning  at  «16  hr,  on  the  average.  The  F-region  scale  height, 
on  the  basis  of  an  idealized  Chapman  profile,  displays  no  consistent  behavior  for  the  three 
months. 

Future  emphasis  will  be  directed  toward  removing  the  seasonal  dependence  from  the 
observations  and  constructing  crossccireiation  functions  between  the  various  ionospheric 
parameters  and  the  magnetic  A  indices  or  their  equivalents.  The  A  indices  are  more  lin¬ 
early  related  to  magnetic  activity  than  are  the  K  indices  which  were  used  in  the  present 
analysis. 
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APPENDIX 

Electron-Density  Distributions  for  March,  April,  and  May  1971 
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Fig.  Al-March  12, 1971 
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Fig.  A2-March  15, 1971 
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Fig.  A3— March  14, 1971 
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Fig.  A6— March  18, 1971 
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Fig.  A7-March  22,  1971 
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Fig.  A8-March  23, 1971 
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Fig.  A9— March  24, 1971 
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Fig.  All— April  21, 1971 
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Fig.  A12-April  22,  1971 


Fig.  AJ3-April  23, 1971 
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Fig.  A14-April  26, 1971 
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Fig.  A15~April  27.  1971 
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Fig.  A16-April  28, 1971 
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Fig.  A17-April  29, 1971 


I 


I 

f{ 

t 

t 

I 

t; 


t 

I 


6 


i 


32  J.  M.  GOODMAN,  M.  W.  LEHMAN,  E.  L.  GOTT,  K.  W.  MORIN,  E.  PIERNIK 


HEIGHT  (km) 


Fig.  AlS-April  30, 1971 
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Fig.  AlO-May  11,  1971 
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Fig.  A20-May  12, 1971 
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Fig.  A21-May  13, 1971 
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Fig.  A22-May  14, 1971 
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Fig.  A23-May  17, 1971 
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Fig.  A24-May  19, 1971 


Fig.  A25— May  20, 1971 
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